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1 Abstract 

I 

I Observational  results  from  shipboard  measurements 

of  temperature  fluctuations,  and  mean  wind,  temperature, 

and  humidity,  are  compared  with  existing  expression  for 

2 

the  surface  flux  and  height  dependence  of  C,p  . Surface 
flux  estimates  are  obtained  from  bulk  aerodynamic  for- 
mula. Temperature  fluctuation  data  are  selected  to 

minimize  a salt-contamination  effect  which  causes  in- 

2 

creases  in  emperature  variance.  Predictions  for  C,^ 
based  on  surface  flux  scaling  agree  within  20%,  except 
for  near  neutral  and  large  unstable  conditions  where  the 
disagreement  can  be  attributed  to  measurement  problems. 


1. 


Introduction 


The  atmospheric  marine  surface  layer  is  that  region  in  the 
first  50  meters  or  so  over  the  sea  which  is  characterized  by 
large  turbulent  vertical  transfer  of  momentum,  heat,  and  mois- 
ture. Progress  has  been  made  in  specifications  of  surface  wind 
fields  and  hydrostatic  stability  for  this  region  since  these 
are  parameters  used  in  estimating  boundary  fluxes  in  numerical 
atmospheric  prediction  and  for  predicting  changes  in  the  upper 
part  of  the  ocean.  Another  interest  in  this  region  is  the 
quantification  of  turbulent  effects  on  optical  wave  propagation. 
The  optical  effects  occur  at  scales  which  are  smaller  than 
those  associated  with  boundary  fluxes.  This  involves  specifica- 
tions of  fluctuations  in  velocity,  temperature,  and  humidity 
for  scales  approaching  those  responsible  for  viscous  dissipa- 
tion of  turbulent  kinetic  energy  and  molecular  diffusion  of 
scalar  (temperature  and  humidity)  quantities. 

Existing  measurement  capabilities  and  turbulence  theory 
enable  observations  of  the  small  scale  properties  of  interest, 
but  the  complexities  and  cost  of  such  direct  observations  in- 
hibit them  on  routine  bases.  Therefore,  there  is  a need  to 
examine  relationships  between  the  small  scale  properties  and 
the  more  feasibly  measured  parameters:  the  mean  winds,  tem- 
perature, and  humidity. 

Descriptions  of  the  small  scale  turbulent  properties  have 
not  been  as  complete  nor  in  the  quantity  for  the  overwater 
regime  as  for  the  overland  regime  (described  by  Champagne,  et 
al.  , 1977),.  The  overwater  regime  requires  separate  considera- 
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I 

! 

tion  because  of  evidence  of  influence  on  the  airflow  by  the  | 

surface  waves  (Davidson,  1974)  and  the  importance  of  the  water 
vapor  content  of  the  air.  The  water  vapor  effects  are  mani-  j 

fested  both  in  atmospheric  stability  and  in  optical  propaga-  i 

tion. 

For  optical  waves,  the  largest  contributors  to  the  refrac- 
tive-index fluctuations  are  temperature  fluctuations  (T').  i 

Friehe,  et  al.  (1975)  have  shown  that  another  significant 
contributor  may  be  the  temperature-humidity  fluctuation  co- 

variance.  They  observed  a contribution  of  24%  to  the  refrac- 

2 

tive  index  structure  function  parameter,  , from  humidity  | 

effects  during  a time  when  the  temperature-humidity  covariance 
was  positive. 

The  intensity  of  temperature  fluctuations  for  scales  from 

a few  millimeters  to  several  meters  can  be  characterized  by 

2 

the  temperature  structure  function  parameter,  C,p  , defined  as 


D,j.(r) 

^2/3 


(1) 


D,j,(r)  is  defined  as  D,j,(r)  = [T(x)  - T(x+r)]^  where  T(x)  and  j 

T(x+r)  are  temperature  at  two  points  on  a line  oriented  normal  j 

to  the  mean  wind  direction  separated  by  the  distance  r.  j 

I 

The  temperature  structure  function  parameter  can  also  be 
related  to  small  scale  processes  of  turbulent  kinetic  energy 
dissipation  and  heat  diffusion  by  the  following  expression 

= n X (2) 


2 


where  n is  a constant  equal  to  3.2,  c is  the  dissipation  rate 
of  turbulent  kinetic  energy  and  x is  the  rate  at  which  the 
temperature  inhomogeneities  are  diffused  or  "smeared"  out. 
Equation  2 arises  from  two  parallel  expressions  for  the  one- 
dimensional variance  spectrum  for  temperature  fluctuations  in 
the  inertial  subrange. 

S.j.(k)  = a k’^/^  (3) 

and 

S.j.(k)  = B X (4) 


where  a and  B are  constants  equal  to  .25  and  .8  respectively. 

The  latter  expression  was  formulated  by  Corrsin  (1951). 

This  paper  describes  results  from  a comparison  of  a very 

2 

substantial  amount  of  C.j,  data  obtained  over  the  open  ocean  with 
theWyngaard  et  al.  (1971)  predictions  based  on  bulk  surface 
layer  measurements.  Friehe  (1976)  performed  a similar  compari- 
son in  which  fewer  data  were  used  and  in  which  a different 

bulk  stability  parameter  was  used.  He  concluded  that  reason- 

2 

able  prediction  of  C,p  could  be  made  from  bulk  formulations. 

2 

2 . Surface  Laver  Bulk  Aerodynamic  Relations  for  C,j. 

2 

Wyngaard,  et  al.  (1971)  related  to  surface  fluxes  with 
Eq.  2 using  separate  expressions  for  e and  x based  on  overland 
measurements.  The  general  expression  was 


z2/3 


f(5):  C 


z 

L 


(5a) 


3 


where 


f(5)  = 


4.9(1  - C < 0 (5b) 

4.9(1  +2.4  5^/^)  C > 0 (5c) 


It  is  noted  that  these  interpolation  formulae  still  lack  sub- 
stantial verification  for  the  overwater  regime. 

These  formulae  agree  with  the  expected  asymptotic  scaling 

2 

(Wyngaard,  1973)  for  under  very  unstable  (local  free 

convection)  and  very  stable  (Z-less)  stratifications.  The 
appropriate  scaling  formulations  in  the  extreme  stability  limits 
are 


free  convection: 


where 


Z-less : 


= constant 


T^  = (TQ^/gZ)^/3  . 


2/3 


constant 


(6) 


w'  9 ' 


(7) 


2 

Therefore,  the  asymptotic  predictions  are  for  C,p  to  vary  as 
-4/3 

Z ' for  free  convection  and  to  be  independent  of  height 
(Z-less)  for  very  stable  conditions. 

In  this  study,  boundary  flux  parameters  in  Eq.  5 , (T^  and 
L)  were  calculated  using  bulk  aerodynamic  formulae.  Bulk 
aerodynamic  formulae  are  those  which  relate  boundary  scaling 
parameters  to  the  wind  speed  at  a level  and  the  temperature 
and  humidity  difference  between  that  level  and  the  surface. 
Their  use  for  indirect  estimates  of  the  boundary  fluxes  has 
received  considerable  attention  in  overwater  observational 
experiments  because  of  difficulties  in  making  eddy  flux,  and 
profile  measurements. 
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The  bulk  aerodynamic  formulae  at  a height  Z for  the  veloc- 


ity and  potential  temperature  scaling  parameters  are: 

u,  - Cp*u  (8) 

T.  = - 63)  (9) 

where  S refers  to  the  surface  value  and  and  Cg  are  drag  co- 
efficients. The  drag  coefficients  depend  upon  height, 
stability,  and  wind  velocity.  Similarly,  for  the  sensible 
heat  flux 

= CjjU  (02  - 0g)  (10) 

where 

C„  - (Cp  Cg)J  (11) 

Since  several  publications  give  detailed  descriptions  of 

the  bulk  method  with  various  experimental  values  for  the  drag 

coefficients,  we  shall  restrict  this  discussion  to  s.  few 

points  that  must  be  considered  when  relating  overwater  measured 
2 

values  of  C,p  to  values  of  T^  and  C based  on  bulk  method 
calculations.  In  particular,  we  must  note  that  T,^  is  directly 
dependent  upon  the  value  of  (Eq.  9).  Using  conventions 
developed  by  Liu  (1978),  the  neutral  stability  temperature  drag 
coefficient  is  given  by  (in  this  convention , (a,j,ic  Z/T^)(39/3Z) 
is  the  dimensionless  profile  function) 

= a,j.K[£n(Z/Z^,j.)]"^  (12) 


r 


i 

! 

i 

[ 

I 

[ 
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where  < is  von  Karman's  constant,  a,p  the  ratio  of  heat  and 
momentum  transport  di f f usivit ies  (at  neutral  stability)  and 
is  the  appropriate  roughness  length  for  the  temperature 
profile.  Since  Wyngaard,  et  al.'s  expression  for  f(^)  was 
based  upon  direct  measurement  of  from  the  heat  flux  and 
Reynold’s  stress,  their  curve  was  not  dependent  upon  a choice 
of  In  view  of  this,  we  will  select  the  value  of  that 

yields  the  best  fit  to  our  data,  treating  this  as  an  indirect 

measurement  of  C-„. 

0N 

The  dependence  of  overwater  neutral  momentum  drag  coeffi- 
cients on  wind  speed  has  been  the  objective  of  many  in- 

vestigations. For  our  analyses,  we  chose  a recently  suggested 
representation  by  Kondo  (1975),  Table  1. 

The  stability  dependence  of  and  Cg  is  defined  on  the 
basis  of  the  following  general  expressions 


= 


'DN 


° [1  - (K)-^ 

c. 


s = 


'0N 


[1  -(a,j,K)"^Cgj,j^  t|'2(C)]^ 


(13) 


(14) 


where  'l'j^(S)  and  stability  corrections  to  wind 

and  scalar  profiles.  Formulations  for  4)^(0  and  4^2^^^  have 
been  presented  by  Businger  (1973)  and  ours  were  similar  except 
for  the  use  of  a,j,  normalization  (Liu,  1978). 

The  dependence  of  the  drag  coefficients  on  stability  leads 
to  computational  considerations  in  estimating  C since  the  co- 
efficients which  define  the  latter  are  themselves  functions 
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of  5 as  indicated  by  the  following  bulk  expression 


(1  - K ^ C^DN 

^ f -1-  -A— ^ 


(1  - (a  ic) 


9N 


where 

Co  i 

^ _ gZ  0N  rA0  + .18AQt 
K = < —!r  r L 5 -* 

° ^ Sn 

All  parameters  are  in  the  MKS  units  with  Q in  gm/kg. 

is  the  equivalent  of  a first  guess  of  ^ based  on  bulk 

differences  and  neutral  values  for  the  drag  coefficients.  For 
each  value  of  wind  velocity,  we  calculated  allowing  cal- 

culation of  C . The  final  value  of  5 was  found  by  solving 
o 

Eq.  15,  iteratively.  An  example  set  of  solutions  is  shown  in 
Fig.  1 for  U = 7.6  ra/sec. 

Friehe  (1976)  effectively  defined  5 as  being  equal  to 

2 

when  he  compared  overwater  C,j,  values  with  bulk  parameter 
predictions  based  on  Eq.  5.  He  used  a constant  value  and 
two  different  values  of  C^j^,  depending  on  the  value  of  UA0  . 

The  differences  between  the  solid  curve  and  dashed  lines  in  Fig. 
1 represent  the  difference  which  would  arise  in  neglecting  stab- 
ility effects  in  the  specification  of  C. 

3.  Experimental  Arrangements 

The  data  were  obtained  during  several  experiments  conducted 
during  a four  year  period  (1974-1978),  over  an  extensive  ocean 
area  and  aboard  three  different  ships  (listed  in  Table  2).  The 
instrumentation  was  designed  for  profile  measurements  of  both 
mean  and  turbulent  parameters.  In  all  experiments,  these  were 


made  at  multi-levels  on  two  masts  separated  spatially  on  the 
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forward  part  of  the  ships  (Fig.  2).  A comprehensive  descrip- 
tion of  the  system  is  the  subject  of  another  paper  (Houlihan, 
et  al. , 1978). 

Mean  temperature  was  measured  with  quartz  oscil- 
lator thermometers,  accurate  to  .01°C.  Mean  humidity  was  mea- 
sured with  Li  Cl  sensors,  accurate  to  about  3IRH. 

Wind  velocity  was  measured  with  cup  anemometers,  accurate  to 

2 

about  5%.  Temperature  fluctuations  (for  C,p  ) were  measured 
with  ac  Wheatstone  bridges  (3.0  KHz  carrier),  using  paired 
resistance  probes  with  2.5um  diameter  platinum  wires.  The 
resistance  probes  were  mounted  on  wind  vanes  to  maintain  correct 
alignment  with  the  relative  wind  direction. 

Two  features  of  the  measurement  arrangement  significant  in 
the  interpretation  of  results  are  the  aspiration  of  mean  tempera- 
ture sensors  and  the  surface  temperature  measurements.  The 
quartz  temperature  sensors  and  humidity  sensors  were  contained 
within  aspirated  weather  shelters.  The  shelters  had  two  purposes 
to  protect  the  sensors  from  the  marine  environment  and  to  reduce 
radiation  effects.  Radiation  effects  were  of  special  concern, 
since  some  of  the  sensors  were  above  the  water  and  some  were 
above  the  deck  of  the  ships,  which  are  heated  by  the  sun  and 
radiate  strongly.  Because  of  radiation  from  the  deck,  we  modi- 
fied the  shelters,  placing  a shield  at  the  base  of  the  unit  to 
reduce  radiation  from  below.  Even  with  these  precautions,  systen 
atic  temperature  corrections  for  the  aspirators  (during  day  and 
night  operations)  have  been  determined  from  calibrations  with  a 
glass  dewar  type  aspirator.  The  daytime  correction  .6C  and  the 
nighttime  correction  is  .4C  (subtracted  from  measured  temperature 
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The  sea  surface  temperature  sensor  (a  quartz  probe) 
was  attached  to  a plastic  hose  so  that  it  floated  within 
one  inch  of  the  surface.  On  all  ships,  its  location  was 
aft  of  the  midship.  The  location  on  all  of  these  ships, 
underway  and  holding  position,  was  definitely  influenced 
by  the  ship's  wake. 

The  mean  data  were  continuously  averaged  and  periodic- 
ally logged  by  a microprocessor  developed  at  the  Naval 
Postgraduate  School  and  adopted  for  the  R/V  ACANIA  by  Plunkett 
(1977).  Mean  values  were  defined  for  15  to  30  minute  periods. 
Most  periods  were  30  minutes  long,  the  planned  duration,  but 
shorter  periods  occurred  due  to  such  events  as  change  in  the 
ship's  course. 

4.  Data  Evaluation  and  Error  Analysis 

2 

Comparison  of  C,p  results  and  available  expressions  are 
made  with  T,^  and  L values  estimated  from  10  meter  wind, 
temperature,  and  humidity  and  surface  temperature  values. 
Initial  plans  were  also  to  estimate  T,^  and  L values  from 
mean  profile  data,  as  evident  from  previously  described 
instru nentation  arrangements.  The  latter  results  are  not 
presented  because  of  inconsistencies  which  have  been  attribu- 
ted to  difficulties  in  defining  small  scalar  and  wind  gradi- 
ents. Although  these  same  errors  existed  in  the  mean  data 
used  for  bulk  estimates,  the  accuracy  requirements  uf  the 
bulk  method  are  much  less  stringent. 
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2 

The  10  meter  mean,  as  well  as  values,  were  obtained  by 
interpolating  from  levels  above  and  below,  if  available,  or  by 
extrapolating  on  the  basis  of  expected  height  variations,  if 
only  one  level  was  available.  With  these  specifications, 
and  C were  computed  for  the  10  meter  level,  using  Eqs.  9 and  15. 

Two  aspects  of  the  measurements  require  further  discussion 
for  this  study.  These  are  the  influence  of  salt  contamination 
on  resistance  wire  sensors,  and  uncertainties  in  the  surface 
temperature  specification.  Although  the  latter  may  not  be  sig- 
nificant for  some  applications,  it  is  significant  in  our  results, 

because  the  air-surface  temperature  difference  occurs  as  a 

2 2 

squared  quantity  in  scaling  C,j,  by  T^  . 

The  influence  of  salt  contamination  of  resistance  wires 

and  thermocouples,  and  hence,  all  temperature  fluctuation 

measurements  in  the  marine  environment,  was  recently  described 

by  Schmitt,  et  al.  (1978).  These  effects  were  discussed  by 

2 

Friehe  (1977)  in  his  comparison  between  overwater  C,j,  values 
and  bulk  aerodynamic  predictions  and  led  him  to  limit  the  data 
considered.  Schmitt,  et  al.  (1978)  suggested  that  erroneous 
temperature  fluctuations  occur  due  to  water  vapor  absorption 
and  evaporation  on  the  salt  nuclei  attached  to  the  wire.  This 
occurs  in  conjunction  with  fluctuations  in  ambient  HH.  The 
effect  is  manifested  by  the  occurrence  of  spikes  in  the  tempera- 
ture trace.  These  spikes  can  cause  overwater  temperature 
variance  spectra  to  exhibit  non-Kolmogorov  slopes. 
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In  an  effort  to  eliminate  the  salt  effect,  the  resistance 

wires  were  washed  or  changed  at  frequent  intervals  during  the 

2 

latest  experiment  (CEWCOM-78).  Mean  ratios  of  C,^,  immediately 

before  and  after  washing  versus  relative  humidity  appear  in 

2 

Fig.  3,  We  observed  that  C,p  values  from  contaminated  wires 
are  larger  for  an  RH  range  from  50%  to  85%,  only.  Presumably, 
the  effect  is  diminished  at  RH  above  85%  because  humidity  fluc- 
tuations are  less  intense  due  to  the  small  RH  gradient.  Their 
effect,  as  described,  is  within  the  system  accuracy.  Presumably, 
the  effect  is  diminished  at  RH  below  50%  because  the  salt  nuclei 
attached  to  the  wires  are  not  activated,  hence,  humidity  fluc- 
tuations do  not  cause  significant  absorption  or  evaporation. 

We  tested  these  interpretations  by  computing  the  ratio  of 
2 2 

measured  C,p  to  computed  C,j,  (using  Eq.  5 and  bulk  parameter 
estimates  of  and  C)  with  data  from  all  experiments  except 
CEWCOM-78.  Mean  values  of  these  ratios  versus  RH  appear  in 
Fig.  4.  The  individual  points  are  mean  values  for  5%  RH 
intervals  and  the  error  bars  are  expected  errors  of  the  aver- 
ages (±  o//W~,  N = number  of  points)  within  each  interval. 

The  number  at  the  top  of  the  error  bar  is  N.  These  data  exhibit 
values  near  1 for  RH  values  above  85%.  There  are  slight  in- 
creases at  95%  and  at  85%  and  there  is  a significant  increase 
at  80%.  The  ratios  are  larger  at  low  RH  values,  compared  to 
corresponding  result  from  clean  wires  in  CEWCOM-78. 

Differences  between  the  results  in  Fig.  3 (washing  and  re- 
placement) and  the  results  in  Fig.  4 (from  the  other  experiments) 
are  the  larger  values  of  the  ratio  at  very  low  RH  values  and  the 

11 


I 


lower  values  above  95%.  The  former  result  is  not  unexpected, 
because  wires  exposed  for  longer  periods  without  cleaning  as 
they  were  for  the  other  experiments  would  be  sufficiently  con- 
taminated to  respond  to  large  humidity  fluctuations,  even  if 
they  were  not  activated.  The  latter  result  could  also  be 
attributed  to  the  degree  of  contamination  and  the  tendency  of 
salt  particles  to  grow  into  water  droplets  at  high  RH.  Con- 
taminated wires  at  large  RH  values  experience  a loss  of  fre- 
quency response,  resulting  in  a lower  sensitivity  to  tempera- 
ture fluctuations. 

Based  on  the  above  results,  we  decided  to  consider  only 

2 2 
those  C,j,  data  measured  at  RH  values  above  82.5%.  C,j,  results 

obtained  with  lower  RH  values  were  excluded  except  from  periods 

immediately  following  sensor  replacement.  All  data  from 

CEWCOM-78  which  were  from  clean  sensors  were  considered. 

2 

Vfe  used  intermediate  C,p  and  A0  results  to  examine  un- 
certainties in  the  sea  surface  temperature.  These  results 

2 

appear  in  Fig.  5 where  C,p  data  were  selected  on  the  basis  of 
RH  considerations  and  A6  values  have  been  adjusted  for  system- 
atic aspirator  corrections.  If  the  measured  surface  tempera- 
tures were  the  "skin"  temperatures,  we  expect  the  minimum  in 
2 

C,p  to  occur  near  A9  equal  to  zero.  This  is  the  case  for  both 

night  and  day.  However,  we  estimate  that  our  A0  values  are 

subject  to  uncertainties  as  large  as  .5°C,  based  on  the  dis- 

2 

placement  and  broadness  of  the  daytime  C,p  minimum.  This  is 
primarily  due  to  an  inability  to  ensure  that  the  sea  surface 
sensor  actually  measured  the  skin  temperature.  Liu  (1978) 
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suggests  how  the  latter  problem  could  have  been  treated  with 
interface  flux  equilibrium  considerations.  However,  his 
suggested  procedures  were  not  applied  to  these  data  because 
they  apply  to  a bucket  temperature  which  was  not  measured  by 
our  sensor. 

The  results  of  an  analysis  of  measurement  errors  appear  in 

Table  3.  These  are  the  errors  introduced  by  instrument  acccura- 

cies,  noise  and  uncertainties  in  certain  calibration  values. 

2 

In  the  case  of  C,j,  , most  of  the  uncertainty  is  due  to  the  salt 
contamination  effect.  In  the  case  of  Tg,  the  sensor  itself 
is  accurate  to  .01°C,  but  one  cannot  be  certain  that  the 
sensor  is  in  fact  measuring  the  surface  temperature. 

5 . Results 

Dimensionless  temperature  structure  parameter  (DTSP)  versus 

C results  appear  in  Fig.  6,  in  the  format  used  by  Wyngaard,  et 

al.  (1971).  Again,  error  bars  are  the  uncertainties  of  the 

means  and  the  niunbers  are  the  number  of  points  defining  the 

mean.  For  this  comparison,  was  determined  by  assuming  that 

Z^,j,  was  independent  of  velocity  and  had  a value  corresponding 

-3 

to  for  a typical  drag  coefficient,  = 1.3  x 10  . The 

2 

appropriate  value  of  for  a measured  C,p  was  estimated  from 
the  following  equation. 


^0N 


= (- 


cl  [1  - (£n  Z/Z^^)-^  t|;,,(5)] 


Trrr 


-) 


A0' 


(17) 


The  value  of  selected  was  that  from  the  CEWCOM-78  data, 
since  these  were  least  influenced  by  salt  contamination. 

Although  it  was  necessary  to  assume  an  initial  value  of 


'0N 


to  obtain  the  values  of  5,  the  final  results  (Fig.  7) 
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i 

j 

I 

are  fairly  insensitive  to  the  initial  choice.  A miniraura  Ic| 
value  of  .05  was  used  in  this  determination  in  view  of  the 

k 

‘ uncertainties  listed  in  Table  3.  The  value  of  selected 

6N 

\ and,  hence,  used  to  compute  T,,,  and  C for  the  results  in  Fig.  6 

1 _3 

\ was  1.3  X 10  . Using  compare  this 

value  with  some  obtained  from  measurements  of  (Table  4),  I 

-3 

assuming  a typical  value  of  1.3  x 10  if  not  given. 

DTSP  results  are  also  presented  with  logarithmic  5 axes, 

Fig.  8.  This  format  enables  examination  of  the  distributions 
of  the  large  number  of  results  near  Z equal  to  zero  in  a more 
illustrative  perspective.  Also,  results  at  ? values  less 
than  -7  which  were  not  included  in  Fig.  6,  can  now  be  viewed. 

This  format  demonstrates  the  exceptional  agreement  between 
the  results  and  the  prediction  over  the  range,  -.03  < 5 < -7. 

Also,  disagreement  between  observed  and  predicted  results  can 
be  readily  related  to  previously  mentioned  measurement  uncer- 
tainties. Disagreements  at  the  minimum  |Cl  intervals,  for 
both  stable  and  unstable  sides,  are  attributed  to  errors  in 
specifying  A0  as  well  as  inherent  noise  in  the  C,p2  measurement 

and  recording  systems.  Disagreement  at  large  |Cl  values  on 

2 

the  unstable  side  is  attributed  to  increased  C,p  values  due 
to  the  salt  contamination  effect,  since  these  data  happen  to 
have  corresponding  relative  humidity  values  near  85%. 

DTSP  results  based  on  asymptotic  scaling  for  unstable 
(Eq.  6,  free  convection)  and  stable  (Eq.  7,  Z-less)  conditions 
appear  in  Figs.  9a  and  9b,  respectively.  Predictions  based  on 
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Eq.  5 (a  and  b)  appear  as  solid  curves, 
the  prediction  was  approximated  as 


For  the  unstable  side, 


z2/3 


9.9 


2/3 


(18) 


which  approaches  a value  of  2.7  for  large  I?].  Eq.  18  is  not 
correct  when  q,^,  as  well  as  , determines  When  q*  is 

negative  (typical  marine  conditions),  the  observed  values 
should  be  less  than  predicted  by  this  expression.  This  is 
the  case  for  these  results.  For  the  stable  side,  the  predic- 
tion is 


C 2 t2/3 

— = 4.9  + 2.4] 


(19) 


which  approaches  a value  of  11.8  for  large  5. 

In  summary,  the  DTSP  results  agree  very  well  with  the 
predictions  under  unstable  conditions  and  reasonably  well 
under  stable  conditions.  Under  stable  conditions,  they  are 
slightly  less  than  the  predictions.  The  determined  value 
agrees  reasonably  well  with  those  reported.  DTSP 
disagreement  under  neutral  conditions  (|C|  < .05)  is  attribu- 
ted to  uncertainties  in  measuring  the  air-sea  temperature 
difference  and  inherent  system  and  measurement  noise.  The 
latter  would  include  the  salt  contamination  effect,  since 
any  humidity  fluctuation  would  cause  apparent  temperature 
fluctuations.  Disagreement  under  very  unstable  conditions, 

C < -7,  was  due  to  the  fact  that  these  data  occurred  with 
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RH  values,  = 85^  , not  totally  immune  to  the  salt  contamina- 
tion effect. 

The  actual  significance  of  the  slight  disagreement  under 
stable  conditions  is  somewhat  uncertain.  However,  it  could 
also  be  attributed  to  the  salt  contamination  effect.  During 
stable  conditions  over  the  ocean,  humidity  and  temperature  j 

fluctuations  are  negatively  correlated.  Hence,  dry-warm  or 
moist-cold  eddies  are  causing  the  fluctuations.  In  the  case 
of  dry-warm  parcels,  contaminated  sensors  would  be  experienc- 
ing a cooling  due  to  evaporation,  but  a compensating  heating 
due  to  the  temperature  change  would  also  occur.  In  the  case 

' 

of  cold-moist  parcels,  the  reverse  would  occur.  j 

DTSP  results  based  on  asymptotic  scaling  appeared  to  j 

approach  the  predicted  constants  at  the  limits.  Agreement  j 

for  these  results  was  expected,  on  the  basis  of  agreement  in  I 

other  DTSP  results,  because  Monin-Obukhov  scaling  was  used  j 

in  determining  T^.  This  arose  when  Q^,  in  Eq.  6,  was  deter-  i 

mined  from  T,^  and  values  which  were  estimated  on  the  basis  ' 

of  Monin-Obukhov  scaling.  i 

2 i 

The  multi-level  C_,  data  enabled  an  examination  of 

^ 1 

predicted  height  variation,  under  different  stability  condi- 
tions. These  results  appear  in  Fig.  10  where  A is  the 

coefficient  obtained  for  the  best  fit  linear  relation  between  i 

I 

2 

observed  Log  C,^,  values  and  Log  Z values. 

Log  = A Log  Z + B.  (20) 
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The  curves  are  predictions  obtained  by  determining  dLog 

2 

C,p  /dLog  Z from  Eq.  5 and  approach  limiting  values  of  -4/3 
and  0 on  the  unstable  and  stable  sides  and  a value  of  -2/3 
at  the  neutral  point. 

Fewer  periods  were  available  for  this  examination  than 
2 

for  the  preceding  results.  This  was  because  probes  at 

upper  levels,  which  were  not  readily  accessible,  were  not 

always  replaced  when  they  became  inoperative.  Also,  we 
2 

measured  C,p  at  a single  level  during  CEWCOM-78  because 

of  frequent  probe  washing  and  replacement. 

Results  in  Fig.  10  indicate  that  height  dependence  of 
2 

observed  C,p  values  agreed  quite  well  with  the  prediction. 

This  is  particularly  the  case  for  the  stability  variation 

of  the  dependence.  A majority  of  the  mean  values  indicate 
2 

that  C,p  decreased  somewhat  less  rapidly  with  height  than 
the  prediction.  This  result  could  be  attributed  to  system 
noise  which  contributed  the  same  at  all  levels  and  would  reduce 
the  apparent  relative  difference  (decrease)  between  levels. 
Evidence  for  this  is  the  larger  (relative)  disagreement  at  near 
neutral  (unstable)  intervals,  where  the  signal  to  noise  ratio 
is  the  smallest. 

6 . Conclusions 

2 

Comparisons  of  a substantial  amount  of  overwater  C,p  data 

have  been  made  with  predictions  of  the  existing  empirical  ex- 

2 

pressions.  Measured  overwater  C,j,  data,  which  were  selected 
to  minimize  the  salt  contamination  influence  are  described 
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very  well  by  Wyngaard,  et  al's  (1971)  expressions,  even 


with  bulk  surface  flux  estimates.  An  estimated  valu 

_3 

, of  1 . 3 X 10  agrees  with  lower  values  obtained  by  others. 

Disagreements  exceeding  measurement  accuracies  (20-30';c)  occur 

2 

for  near  neutral  conditions  where  is  insignificant  in  most 

applications,  and  for  large  unstable  conditions.  These  dis- 
agreements are  attributed  to  measurement  difficulties  rather 
than  to  the  validity  of  the  existing  expressions.  Agreement 
between  observed  and  predicted  height  distributions  over  the 
full  stability  range  provide  further  verification  of  the  e.xist- 
ing  expression. 

In  conclusion,  a typical  (but  not  the  best)  comparison  of 

2 

observed  and  predicted  C,j,  values  during  a 24  hour  period  in 
CEWCOM-78  is  presented  in  Fig.  11a  to  further  demonstrate  the 
overwater  applicability  of  the  e.xpressions . The  predicted 

values  based  on  bulk  estimates  of  surface  fluxes  are  within  i 

the  accuracy  of  the  measurements  in  describing  the  observed 
I values.  Coincident  bulk  |5|  estimates  for  the  period  appear 

^ in  Fig.  lib. 

I 
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Table  1, 


^DN  wind  speed  (10  meter) 

from  Kondo  (1975). 


U ms 


-1 


.3  - 2.2 
2.2  - 5.0 

5.0  - 8.0 

8.0  - 25.0 


1.08  X U 
.77  + .086  X U 
.87  + ,067  X U 
1.2  + .025  X U 


TABLE  2.  Summary  of  Experiments 


NUMBER  OF 


EXPERIMENT 

SHIP 

LOCATION 

PERIOD 

PERIODS 

R/V  ACANIA'* 

Monterey  Bay 

All  Months 
1973-76 

98 

CEWCOM-76^ 

R/V  ACANIA 

E.  Pacific 

Sep-Oct  76 

285 

NRL  - EOMET^ 

USNS  HAYES^ 

N.  Atlantic 
& Mediterranean 

May-Jun  77 

269 

USNS  KANE® 

Mid-Atlantic 

Feb-Mar  78 

38 

CEW COM-78^ 

R/V  ACANIA 

E.  Pacific 

May  1978 

200 

Total 

890 

^Cooperative  Experiment  West  Coast  Oceanography  & Meteorology-1976 
^Naval  Research  Laboratories,  Electro-Optical  Meteorology  Cruise 
^Cooperative  Experiment  West  Coast  Oceanography  & Meteorology-1978 

4 

Operated  by  Department  of  Oceanography,  Naval  Postgraduate  School, 
Monterey,  California 

5 

Operated  by  Naval  Research  Laboratories,  Washington,  D.  C. 

Operated  by  Naval  Oceanographic  Office,  Washington,  D.  C, 
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Table  3.  Error  analysis  for  measured , 

intermediate , and  scaling  parameters 


Measured 


T T 

S 


U 


RH 


RMS  Error  (±)  . 3°C  . 2°C  8% 


3% 


30% 


Intermediate 


Ae 


AQ 


Accuracy  (1) 


.5°C 


.4  gm/Kg 


Table  4. 


r. 


Empirically  determined  sensible 
heat  drag  coefficients, 

X lO’ 


Sources 

c 

HN 

c 

^DN 

^6N 

Pond,  et  al.  (1974) 

1.5 

1.45 

1.55 

Dunckel , et  al . (1974) 

1.5 

1.60 

1.40 

Muller-Glewe  and  Hinzpeter  (1974) 

1.0 

(1.30)^ 

.77 

Smith  (1974) 

1.2 

1.02 

1.41 

Friehe  and  Schmitt  (1976) 

.91 

^ (1.30) 

.64 

This  study 

1.  3 

1.30 

1.  34 

) denotes  assumed  Cj^  value;  other  given  by  authors 

2 

Excludes  high  wind  data 
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Figure  1. 


Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 


Ratio  of  C based  on  stability  dependent  drag 
coefficients  (5)  to  C based  on  constant  drag 
coefficients  (C  ) versus  C . 


2.  Schematics  of  shipboard  mounting  arrangements. 

(a)  R/V  ACAx\IA,  (b)  USNS  HAYES,  and 
(c)  USNS  KANE. 

2 

3.  Ratio  of  C,p  (before  to  after  replacement  or 
washing)  versus  RH;  (From  CEWCOM-78). 

2 2 

4.  Ratio  of  measured  Cj,  to  computed  C.p  (from  Eq.  5) 
versus  RH,  for  all  except  CEWCOM-78  data. 

5.  versus  A9  (a)  daytime  (0800-1900  LST), 

(b)  nightime  (1900-0800  LST). 

6.  Mean  DTSP  versus  C results.  Solid  curve  is 
prediction  based  on  Eq.  5. 

7.  Cgj^ versus  C results.  C^^^  based  on  Eq.  17  and 
CEWCOM-78  results. 

8.  Same  as  Figure  6 with  logarithmic  5 axis. 

9.  DTSP  based  on  asymptotic  scaling  (a)  free  convection 
curve  is  approximate  prediction,  Eq.  18  and, 

(b)  Z-less,  curve  is  prediction,  Eq.  19. 

2 2 

10.  Height  variations  of  C,j,  versus  C.  A(C,j,  ) is  defined 
by  Eq.  20,  solid  curves  are  prediction  based  on  Eq.  5. 

11.  Results  from  24  hour  period  during  CEWCOM-78 

(a)  Observed  and  Predicted  C^  for  10  meters 

(b)  Observed  C- 
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Figure  2, 


Schematics  of  shipboard  mounting  arrangements, 
(a)  R/V  ACANIA,  (b)  USNS  HAYES,  and 
(c)  USNS  KANE. 
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Ratio  of  (before  to  after  replacement  or 

washing)  versus  RH;  (From  CEWCOM-78). 


H ro 


Figure  9.  DTSP  based  on  asymptotic  Scaling  (a)  tree  cc 
curve  is  approximate  prediction,  Eq . 18  and, 
(b)  Z-less,  curve  is  prediction,  Eq . 19. 
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Figure  11.  Results  from  24  hour  period  during  CEWCOM-78 

(a)  Observed  and  Predicted  C,p  for  10  meters 

(b)  Observed  C. 
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